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The line-of-sight distribution of water in the SgrB2 complex 
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Abstract. We report the detection, with the Caltech Submillimeter Observatory*, of the 894-GHz HDO(li,i — 0o,o) 
transition, observed in absorption against the background continuum emission of the SgrB2 cores M and N. 
Radiative transfer modeling of this feature, together with the published data set of mm and submm HDO and 
H3> 8 transitions, suggests that ground-state absorption features from deuterated and non-deuterated water trace 
different gas components along the line of sight. In particular, while the HDO line seems to be produced by the 
large column densities of gas located in the SgrB2 warm envelope, the H3, 8 ground-state transition detected by 
SWAS and KAO at 548 GHz (Neufeld et al. 120001 Zmuidzinas et al. ll995al is instead a product of the hot, diffuse, 
thin gas layer lying in the foreground of the SgrB2 complex. 
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1. Introduction 

The SgrB2 cloud is one of the most massive star form- 
ing regions in our Galaxy. It has many unique characteris- 
tics, among them an exceptional chemistry. Several species 
(FeO, Walmsley et al . 1201121 NH 2 , van Dishoeck et al. lT5551 
HF, Neufeld et al. H5571 NH, Goicoechea, Cernicharo & 
Caux 2000) have, in spite of searches elsewhere, been de- 
tected only toward this source. One possible explanation 
for this enigmatic chemistry is the existence of a fore- 
ground layer of hot gas, which has been observed e.g. in 
high-excitation ammonia transitions (cf. Huttemeister et 
al. 11995(1 . However, since this layer has the same velocity 
as the ambient gas in SgrB2, it has been proven difficult 
to assess its importance for the chemistry of many species. 
In particular its importance for the water chemistry has 
been a matter of recent debate (Ceccarelli et al. 2002 , here- 
after C02; Neufeld et al., |in prep~| ). While its temperature 
and density are well constrained by observational data (cf. 
Huttemeister et al. I1995fl . its column density and spatial 
extent have remained elusive. In this paper we are able to 
determine its importance for water chemistry and provide 
constraints on its column density and spatial width, by 
modeling the HDO and H^O emission and absorption. 

Water is known to be a fundamental ingredient of 
the interstellar medium. It is a major coolant of star- 
forming clouds (Ceccarelli, Hollenbach & Tielens, ^996) 
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and therefore it affects the dynamical evolution of the 
clouds. A ubiquitous tracer of shock-heated gas, it dra- 
matically influences the chemistry in shocked regions (cf. 
Neufeld & Melnick lTWl Bergin, Melnick & Neufeld lT55g|l . 
Although direct, ground-based observations of non-masing 
water lines are made extremely difficult by atmospheric 
absorption, H2O abundances can be estimated via obser- 
vations of isotopomers, such as HJ, 8 (e.g. Phillips et al. 
1T5781 Jacq et al. 115551 and IT5501 Gensheimer et al. lT55Bjl . 
and, in recent years, by the availibility of satellites such 
as ISO (e.g. van Dishoeck & Helmich 1996; Cernicharo 
et al. 1T5571 Wright et al. I2000|l and not ably SW AS (cf. 
Melni ck et a l. 120001 Snell et al. l2()()()al and l2000bl Neufeld 
et al. 120001 hereafter N00), but also by airborne obser- 
vations of H-^O (e.g. Zmuidzinas et al. 11995a! hereafter 
Z95a; Timmermann et al. 1996). However, the deuterated 
counterpart of water, HDO, presents many features ob- 
servable from ground in the cm, mm and submm wave- 
length atmospheric windows (cf. Henkel et al. 119871 Jacq 
et al. 1T550I Schulz et al. 1T55T1 Helmich, van Dishoeck & 
Jansen lT556l Jacq et al. H5551 Pardo et al.|200TJ), and can 
also be used as a tracer of water abundance under the 
assumption that both the deuterated and non-deuterated 
species are spatially coexistent, and that their abundance 
ratio is constant throughout the region of interest. 

We report here the detection, with the Caltech 
Submillimeter Observatory, of the ground-state (li,i — 
0o,o) transition of deuterated water, observed in absorp- 
tion against the background continuum sources SgrB2(M) 
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and SgrB2(N). The radial velocity of the observed HDO 
feature suggests their direct connection to the SgrB2 com- 
plex. 

In § 13.11 we estimate the column density of the ab- 
sorbing HDO using classical absorption line assumptions, 
based on the 894-GHz feature. Moreover, thanks to the 
H^ 8 SWAS data of N00, we are able to estimate the 
[HDO]/[H 2 0] ratio towards the SgrB2 complex (§ 
The interpretation of these results, however, depends crit- 
ically on the assumptions on the location of the absorbing 
gas. In fact, the observed absorption could be produced 
both in the above-mentioned hot gas layer, and in the 
warm envelope of molecular gas in which the two main 
cores are known to be embedded. Moreover, the assump- 
tion that deuterated and non-deuterated water be spa- 
tially coexistent is probably incorrect in the SgrB2 cloud. 
In § ^] we make use of our CSO 894-GHz data, as well as 
of the published HDO and H 2 8 observations carried out 
at mm and submm wavelengths, to model the distribution 
of water throughout the whole SgrB2 complex in its three 
component: hot cores, warm envelope and hot layer, using 
a state-of-the-art radiative transfer model. The results are 
discussed in § 



2. Observations 

The observations were carried out at the Caltech 
Submillimeter Observatory atop Mauna Kea, Hawaii, 
on July 28, 2001. Scans were taken, using the chop- 
ping secondary with a throw of 4', towards the coor- 
dinates a J20 oo = 17 h 47 ro 20?206, <5 J2 ooo = -28°23'05'/27 
for SgrB2(M), and a J2 Q00 = 17 h 47 m 20!389, 6 J20Q0 = 
-28°22'22'.'25 for SgrB2(N). At the observing frequency 
of 893.6 GHz, the CSO 10.4-m antenna has a HPBW of 
about 10". 

The 500-MHz facility AOS was used as backend, pro- 
viding a velocity resolution of ~ 0.5 kms" 1 . System tem- 
peratures varied between 7500 and 8000 K. Several dif- 
ferent LO settings were needed in order to cover a wide 
enough velocity range, and to access at least two, namely 
the 65- and 81-kms _1 , of the several velocity components 
known to exist in the gas clouds lying on the line of sight in 
the direction of the two cores (e.g., Whiteoak & Gardner 
IW9l Martin-Pintado et al. UMl Greaves et al. IT532I 
Tieftrunk et al. [1994). Also, the use of various LO set- 
tings allowed us to rule out possible contamination from 
signal coming from the image sideband. In total, the range 
-eOkms^ 1 <> 

vlsr ^ 140 km s was covered. However, 
the low signal-to-noise ratio does not allow an analysis 
of the other velocity components, e.g. at ~ and ~ 100 
km s _1 , towards which water absorption has been detected 
(cf. N00). 

We would also like to mention that the second ground- 
state rotational transition of HDO, the lo,i-0o,o at 465 
GHz, has been searched for by us as well as by other au- 
thors (e.g., E. Bergin and collaborators, priv. comm.). A 
clear detection could not be achieved, however all data 



sets are consistent in indicating that the 465-GHz transi- 
tion shows up as a weak emission line. 

Pointing is a delicate issue for the 850-GHz receiver 
(Kooi et al. 2000) at the CSO, therefore it was checked 
rather often through 5-point maps of Mars, located only 
~ 11° away from our target sources. The accuracy of our 
pointing is confirmed by the continuum levels we measure 
towards the two cores (21 K for M and 13 K for N, in 
units of main-beam temperature), that match, within 20% 
error, with the values inferred by a map of the 350- /im 
continuum emission acquired with SHARC 1 at the CSO 
(Dowell et al. I1999f) . Observing Mars also allowed us to 
measure the beam efficiency of the telescope, which we 
found to be around 30%. 

3. Results 

3.1. Determination of HDO column densities 

Fig- El displays two spectra of the ground-state li,i — 0o,o 
transition of HDO, observed in absorption against the con- 
tinuum emission of SgrB2(N) (upper panel) and SgrB2(M) 
(lower panel). 

The determination of column densities from absorption 
lines is in general more accurate, with respect to the quan- 
tities inferred from the analysis of emission lines, because 
the calculated value does not depend on the excitation 
temperature of the line, T ex , as long as this is negligible 
with respect to the temperature of the background con- 
tinuum source, Tc: in other words, no assumptions need 
be made about the physical state of the gas, except that 
T cx <C Tc and that only the ground state of the molecule 
is populated (we will show in § 0] that neither holds for a 
large part of the warm envelope). The total column den- 
sity of the absorbing species is given by: 



AT 87W a 91 



(1) 



The opacity of an absorption line can be determined by: 
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Tl 



(2) 



where Tl/Xc is the line-to-continuum ratio. Therefore, 
eq. Q only depends on T^/Tc and on the line width, Av, 
quantities that can be easily derived from a Gaussian fit of 
the optical depth. Note that, at this stage, no assumptions 
have been made about the actual location (warm envelope 
or hot layer) of the gas producing the observed features. 

The fit to the 65- and 81-kms^ 1 components of N, and 
to the 65-kms _1 component of M are shown in Fig. ^ 
while the derived parameters are summarized in Tab. ^ 
The relative values of the total column density are listed 
in Tab. [21 Studies of the 1.3-mm dust emission in the 
region (IRAM 30-m maps, Gordon et al. [1993) indicate 
that the column density of H 2 , iV(H 2 ), reaches values of 
2.2 x 10 24 cm" 2 towards SgrB2(M), and 7.5 x 10 24 cnT 2 
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towards SgrB2(N). Such results agree, within a factor of 2, 
with previous estimates obtained by Goldsmith, Snell & 
Lis (1987) based on lower-spatial-resolution 1.3-mm con- 
tinuum maps. The relative abundance of gas-phase HDO 
displayed in the last column of Tab.|21has been estimated 
on the basis of the above mentioned values of 7V(H2). The 
estimate takes into account that the amount of gas we 
"see" through the HDO absorption is only a half of the 
actual amount of gas measured through the core. We take 
the error on the HDO abundance to be of ~ 50%. 
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Fig. 1. The 894-GHz HDO (1 M - ,o) transition, ob- 
served in absorption against the continuum emission of 
SgrB2(M) (lower panel) and of SgrB2(N) (upper panel). 
The y-scale gives the line-to-continuum (T^/Tc) ratio. 
The absorption towards SgrB2(M) shows one velocity 
component only, around 65 kms -1 , while the gas in front 
of SgrB2(N) shows an additional component around 81 
kms -1 . Such behaviour has been found for other molecu- 
lar species as well (e.g. H2CO, see Martin-Pintado et al. 
1990). The Gaussian fits to the optical depth of the lines 
are shown as solid curves, and the parameters derived from 
the fit are summarized in Tab. 



Table 1. Results of the fits to the optical depth of the 
HDO absorption features shown in Fig. ^ Columns 2 to 
4 give, respectively, the radial velocity of the absorbing 
gas, the line opacity and the line width. The fits have 
been performed with the method ABSORPTION from the 
GILDAS CLASS package. 



SOURCE 


VLSR 


T 


Av 




(kms" 1 ) 




(kms -1 ) 


SgrB2(M) 


64.7 ±0.5 


0.38 ±0.04 


12.5 ± 1.3 


SgrB2(N) 


63.7 ±0.5 


0.77 ±0.05 


14.6 ± 1.3 


SgrB2(N) 


82.1 ± 0.9 


0.32 ± 0.04 


13.0 ±1.7 



Table 2. Column densities of HDO calculated from the 
absorption features shown in Fig. ^ Column 2 gives the 
radial velocity of the absorbing gas component. The values 
of the column density (column 3) are derived, using eq. (JJJ, 
assuming that only the ground-state level of HDO is pop- 
ulated. Column 4 lists the estimated HDO abundance, cal- 
culated assuming a H2 column density of 2.2 x 10 24 cm -2 
towards SgrB2(M) and 7.5 x 10 24 cm -2 towards SgrB2(N) 
(Gordon et al. ITM^ . The estimated error of the HDO 
abundance is of order 50%. 



SOURCE 


VLSR 

(kms -1 ) 


N(HDO) 

(xlO 13 cm -2 ) 


[HDO] 
(xlO" 11 ) 


SgrB2(M) 


64.7 


1.2 ±0.2 


1.1 


SgrB2(N) 


63.7 


2.8 ±0.3 


0.7 


SgrB2(N) 


82.1 


1.1 ±0.2 


0.3 



3.2. The [HDO]/[H 2 0] ratio 

N00 have observed, with the Submillimeter Wave 
Astronomy Satellite (SWAS), the ground-state lj. — lo.i 
transition of o-H^O (547.7 GHz), in absorption against 
the continuum emission of SgrB2. They were therefore 
able to estimate the column density of H^O over a wide 
range of radial velocities (from -120 to 20 kms -1 ). 

We use the H^O spectrum of N00 to estimate the 
[HDO]/[H2 0] ratio in the gas components with radial ve- 
locities of 63 and 81 kms -1 . The comparison between the 
two data sets is not straightforward, since the spatial res- 
olution of SWAS, 3f 3 x 4.'5, is such that the emission from 
SgrB2(M) and SgrB2(N), which are ~ 47" apart, is not 
resolved. To get around this problem, we simply sum up 
the HDO scans towards M and N weighted by the atten- 
uation of the circularized 3f 9 SWAS beam, and compare 
the resulting spectrum (hereafter M±N) to that of N00. 
Such procedure is based on the assumption that the ab- 
sorbing cloud is extended and uniform over the SWAS 
beam size. This is a reasonable approximation for the 65- 
kms -1 component, which is known to be extended over 
both SgrB2(M) and N. However, the 81-kms -1 gas com- 
ponent is definitely smaller than the SWAS beam (it is 
only observed towards the northern core). 

Fig. |U compares the HDO absorption features towards 
M+N (upper panel) with the 548-GHz o-H^O ground- 
state line observed by N00 (lower panel) . The solid curves 
in Fig. [2 represent the two-component Gaussian fits to 
the opacity of the observed features, and the fitting pa- 
rameters are summarized in Tab. [3] Again, we used eq. £Q| 
to calculate the total column density of the two species, 
assuming that, for both of them, only the ground-state 
level is populated. The results are summarized in Tab. 
The last column of Tab. 21 shows the estimated value of 
the [HDO]/[H20] ratio for the two velocity components of 
the gas, calculated assuming that [ 16 0]/[ 18 0] = 261 ± 20 
(Whiteoak & Gardner HHED>, and that [o-H 2 0]/[p-H 2 0] = 
3. We find the [HDO]/[H 2 0] ratio towards SgrB2, at 63 
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and 81 kms -1 , to be ~ 5 x 10~ 4 and ~ 10~ 3 respec- 
tively (see Tab.rjJ). These values are about 30 and 70 times 
higher than the measured deuterium abundance in the 
Local Interstellar Medium (1.5xl0~ 5 , Linskv lT99gjl . The 
enhancement of the deuterium fractionation observed in 
the [HDO]/[H20] ratio appears even more important if 
one considers that deuterium has been measured to be 
under-abundant in the Galactic Center region. Lubowich 
et al. ( 2000) have estimated the [D]/[H] ratio towards SgrA 
to be around 1.7 x 10~ 6 . Similar values have been measured 
towards SgrB2 by Jacq et al. (1999 ) and by Polchampton 
et al. ( 2002 ) . All estimates are affected by very large uncer- 
tainties, however an order-of-magnitude comparison with 
our measured [HDO]/[H 2 0] ratio shows that it is actually 
a few 10 2 times higher than the deuterium abundance in 
the region. This value, scaled to the lower [D]/[H] ratio 
in the Galactic Center, is consistent with that predicted 
by the steady-state chemical models of Roberts & Millar 
(2000) for gas temperatures ranging between ~ 30 and 
100 K and densities between ~ 10 3 and 10 s cm~ 3 , roughly 
the range of temperatures and densities expected in the 
warm envelope. However, the analysis presented in the 
next section shows that the HDO and H2O absorption are 
actually produced in different locations, hence the above 
mentioned value of the [HDO]/[H20] ratio has little rele- 
vance. We will address the issue again in § 

4. Modeling the distribution of water 

As previously mentioned (§ Q, it is known from obser- 
vations of, e.g., H2CO (cf. Whiteoak & Gardner 119791 
Martin-Pintado et al. 1990) that the two main cores in the 
SgrB2 complex - SgrB2(N) and SgrB2(M) - are embedded 
in a warm (15 K ^ T gas < 100 K), dense (n(H 2 ) ~ 10 5 
cm -3 ) envelope of molecular gas. Also, a hot (500 K <, 
T gas ^ 700K) and diffuse (10 3 c m - 3 < n(H 2 ) £ 10 4 c m - 3 ) 
gas layer is located in the foreground, as indicated by the 
detection in absorption of high-energy ammonia transi- 
tions (Wilson et al. 119521 Hiittemeister et al. H99B1 C02). 
Both gas components show characteristic radial velocities 
around 65 and 81 km s -1 (the latter being observable only 
towards the northern core), hence the hot foreground gas 
layer is thought to be physically connected to the cloud 
complex as well. 

The interpretation of the results illustrated in § 13.21 
strongly depends on the assumptions about the location 
of the absorbing gas. Absorption due to molecules such 
as H2O and its isotopomers has generally been attributed 
to the warm envelope (Z95a; N00), under the assumption 
that only the ground-state level was populated. However, 
C02 have proposed that all the absorbing water might in- 
stead be located in the hot gas layer. In this case, accord- 
ing to C02, the estimated column density of H2O towards 
the main velocity components would be about one order 
of magnitude higher than that quoted, for example, by 
Z95a and by us in § 13.21 Such a high value would support 
the C-shock model used by Flower and collaborators to 
explain the observed features of the hot layer, such as the 
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Fig. 2. The 894-GHz HDO(l M - 0)0 ) and the 548-GHz 
o-H2 8 0(li i o — lo.i) absorption features observed towards 
SgrB2. The spectrum in the upper panel is the sum of the 
scans relative to the HDO absorption towards SgrB2(M) 
and SgrB2(N) separately (see Fig. observed at the 
CSO with a 10" beam. The lower panel shows the o-H^O 
absorption towards SgrB2, observed by N00 with SWAS 
(3'3 x 4.'5 beam). The y-scale gives the line-to-continuum 
(Tl/Tc) ratio. The solid curves represent the Gaussian 
fits to the line opacity. The radial velocities, line opaci- 
ties and widths, derived from the fit for both transitions, 
are listed in Tab. El The larger width of the SWAS o- 
H5, 8 line, with respect to the HDO feature resulting from 
the combination of the CSO scans, agrees with the ex- 
tended water distribution observed by, e.g., Cernicharo et 
al. (1997) and Neufeld et al. ( |in prep.p . 

high gas temperature, low gas density and limited spatial 
width (Flower & Pineau des Forets 1T994I Flower, Pineau 
des Forets & Walmslev 1T995) . 

From a chemical point of view, it is entirely possible 
that a large percentage of the observed water column den- 
sity be located within the diffuse hot layer rather than in 
the warm envelope. In fact, the gas temperature in the hot 
layer is sufficient not only to induce evaporation of water 
ices from the surface of dust grains, but also to trigger 
the gas-phase production of H2O via the following chain 
of neutral-neutral reactions: 

+ H 2 -» OH + H, 
OH + H2 -» H 2 + H, 

which proceed very rapidly as T <; 400 K (Elitzur & 
de Jong lT978l also cf. Bergin, Melnick & Neufeld [T99j§ . 
Therefore, we can expect the abundance of H2O to be 
enhanced in the hot layer with respect to the warm enve- 
lope, up to the point where all the oxygen not locked in 
CO is converted into H2O. In fact, HDO can be formed in 
shocked gas as well, via the reactions (Bergin, Neufeld & 
Melnick (JT9J9J: 

O + HD -> OD + H, 
OD + H2 -> HDO + H, 
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Table 3. Results of the Gaussian fits to the optical depth of the HDO (columns 1-3) and o-HJ, 8 (columns 4-6, N00 
data) absorption features observed towards SgrB2 and shown in Fig. |2 The listed parameters are, for each transition, 
respectively the radial velocity, the line opacity and the line width. The fits have been performed with the method 
ABSORPTION from the GILDAS CLASS package. 



hdo(i m - o , ; 


) 


o- 


Hf O(l ll0 - 1 


0,l) 


VLSR 


T 


Av 


VLSR 


T 


Av 


(kms" 1 ) 




(kms" 1 ) 


(kms -1 ) 




(kms -1 ) 


64.3 ±0.3 


0.49 ± 0.03 


13.3 ±0.7 


61.9 ±0.3 


0.77 ±0.02 


17.7 ±0.7 


82.0 ±0.0 


0.11 ± 0.02 


14.6 ± 0.2 


80.3 ±0.9 


0.13 ±0.02 


9.6 ±2.7 



Table 4. Total column densities of HDO and o-H 2 8 in the 3.'3 x 4.'5 SWAS beam, producing the absorption features 
shown in Fig. [2 Column 1: average radial velocity of the gas component (calculated from the values in Tab. 01 columns 
1 and 4); column 2 and 3: total column density of HDO and o-H^O calculated using eq. 2] and the fit parameters in 
Tab. |3 Since we assumed only the ground-state level the absorbing material to be populated, these are to be considered 
as lower limits. Column 4: [HDO]/[H 2 0] ratio, calculated on the basis of the estimated column densities of HDO and 
o-H^ 8 assuming that [ 16 0]/[ 18 0] = 261 ± 20 (Whiteoak & Gardner EHETJ) and that [o-H 2 0]/[p-H 2 0] = 3. 



vlsr N(HDO) 
(kms" 1 ) (10 13 cm" 2 ) 


N(o-HfO) 
(10 13 cm" 2 ) 


N(H 2 0) 
(10 16 cm" 2 ) 


[HDO]/[H 2 0] 

(io- 3 ) 


~ 63 1.67 ±0.13 
~81 0.40 ±0.07 


6.33 ± 0.30 
0.60 ±0.19 


2.20 ±0.20 
0.21 ±0.07 


0.8 ±0.1 
1.9 ±0.7 



and 

OH ± HD -> HDO ± H. 

But, because of the lower HD abundance, the timescales 
are slower for these reactions than for the H 2 pro- 
duction paths, such that all the available atomic oxy- 
gen will be locked in H 2 prior to the formation of 
HDO. Thus, gas-phase chemistry in the hot layer yields 
[HDO]/[H 2 0] ~ [D]/[H]. We believe this picture cannot 
be significantly altered by grain chemistry, because i) the 
temperature of the dust in the hot layer is known to be 
much lower than that of the gas (cf. Wilson et al. 1982 
Huttemeister et al. I1993|) . too low to allow evaporation of 
ice mantles from the grain surface; ii) sputtering of dust 
grains, on the other hand, cannot account for the release 
of intact molecules; and however, Hi), the build-up of ice 
mantles is not efficient at low densities, hence the deposit 
of water ices on the grains surface is likely to be negligible 
in the first place 2 . 

To summarize, we believe that it is reasonable to ex- 
pect the contribution of the hot layer to the observed HDO 
absorption to be negligible, whereas its contribution to the 
observed H 2 absorption is likely to be significant. It is 
important to estimate the extent of such a contribution: 
a precise determination of the column density of water 

2 This is supported, for example, by observations of high 
abundances of gas-phase atomic oxygen in molecular clouds 
with density ~ 10 2 - 10 3 cm -3 (Lis et al. 12001^ . suggesting 
that much of the oxygen which is not locked in CO could well 
be present in the gas in its atomic form, rather than depleted 
on dust grains. 



would set a tight constraint on the physical models that 
attempt to identify the heating mechanism responsible 
for the high temperatures in the hot layer. Unfortunately, 
since both gas components show roughly the same radial 
velocity, it is not possible to separate, by purely observa- 
tional means, the ground-state water absorption produced 
in the warm envelope from that produced in the hot layer. 
However, a number of HDO and H^O transitions have 
been observed towards the SgrB2 cores with a variety of 
instruments (see Tab.|S]and Fig.[3J), and detailed modeling 
can be performed to disentangle the contribution to the 
observed features from the different cloud components. 

The available dataset, enriched by our detection of 
the ground-state HDO transition at 894 GHz, provides 
sufficient observational constraints to model, in a self- 
consistent manner, the HDO and H 2 abundance in all 
three components of the SgrB2 cloud. In detail: 

a) the observed emission lines of p-H 2 8 at 203 GHz 
(SgrB2(N), Gensheimer et al. iL^jl . and of HDO at 
143, 226, 241 GHz (SgrB2(M) and N, Jacq et al. rrWIjl 
constrain the HDO and water abundance and the H 2 
density in the hot-core-type components; 

b) being virtually absent in the hot layer, HDO can be 
used to calibrate the water content of the warm enve- 
lope. In particular, the 894-GHz HDO absorption fea- 
ture (this work) observed towards both SgrB2(M) and 
N, together with the [HDO]/[H 2 0] ratio estimated at 
point a), help us set the HDO and H 2 abundance in 
this region; 

c) finally, having estimated the column density of water 
in the warm envelope, it is possible to assess whether or 
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Fig. 3. Rotational levels of deuterated water up to ~ 
400 K. The transitions included in our data set are in- 
dicated (solid lines: this work; dotted lines: Jacq et al. 
IT35U1 see Tab.0. 

not the quantity of non-deuterated water in the warm 
envelope can be responsible for the HJ, 8 absorption 
observed at 548 GHz (N00, Z95a). The 18 0-^m o-H^O 
line observed by Cernicharo et al. I|1997p towards the 
SgrB2 complex is blended with the H30 + (17/-lf ) line 
(Goicoechea & Cernicharo 2001), and cannot therefore 
be used as a reference for our purposes. 

4.1. The model 

We use the static radiative transfer code described by 
Zmuidzinas et al. I1995bl (hereafter Z95b) to reproduce the 
intensities observed for the features listed in Tab. [5] A few 
changes have been made with respect to the original ver- 
sion, but the bulk of the model is the same and can be 
summarized as follows: 

— A background, internally heated continuum source is 
embedded in a spherical molecular envelope of radius 
0.05 pc< r < 22.5 pc (corresponding to the warm en- 
velope as defined at the beginning of §0J. 



— The large-scale structure of the cloud is well con- 
strained by the observations of C 18 carried out by Lis 
& Goldsmith iJIMfl hereafter LG89). The density pro- 
file adopted in the model, as well as the radial variation 
of the dust temperature, correspond to their Model C. 
The temperature of the gas is coupled to that of the 
dust, Tdust = T gas = T. 

— The abundance of the modeled species, assumed to be 
constant throughout the envelope, is a free parameter, 
as well as the width (in kms -1 ) of the modeled feature. 

— The code allows radiative transfer calculations for one 
Gaussian velocity component only. In what follows, we 
will only try to reproduce the 65-kms _1 component of 
the molecular gas. 

This geometrically straightforward model has proved 
successful in reproducing the dust continuum emission 
observed towards the SgrB2 cloud complex at mm and 
submm wavelengths, as well as absorption features such 
as the HCl(l-O) and H^Ofli^o — lo.i) lines observed in 
the submm band with the Kuiper Airborne Observatory 
(KAO, see Z95a and Z95b). Nevertheless, it is not capable, 
as it is, to reproduce those HDO and H^O features that 
have been observed in emission, and that are thought to 
be produced in the hot-core- type components of SgrB2(M) 
and SgrB2(N) (cf. Jacq et al. 119901 Gensheimer et al. 
1996). In fact, a molecular hot-core component, intended 
as a compact, dense, warm molecular region character- 
ized by fairly flat radial profiles of density and temper- 
ature, is missing in the model, since the region of ra- 
dius r < 0.05 pc is devoid of molecular gas. A radius 
of ~ 0.05 pc has indeed been attributed to the hot cores 
embedded in SgrB2(M) and N, based on the high-spatial- 
resolution interferometric maps of the continuum emission 
towards the complex (cf. Lis et al. 119931 hereafter L93). 
Therefore, we have modified the model by introducing a 
region of radius r < 0.05 pc, in which the H2 density and 
the gas and dust temperatures are constant and equal 
the values of density and temperature at r = 0.05 pc, 
based on the radial profiles of LG89, Model C. This yields 
n(H 2 ) ~ 3.4 x 10 7 cm" 3 and T dust = T gas ~ 200 K in 
the hot-core component. Both values are in good agree- 
ment with the estimates based on the continuum emission 
(n(H 2 ) ~ 3 x 10 7 cur 3 in the hot core of SgrB2(M), and 
n(H 2 ) - 2 x 10 7 cm~ 3 in SgrB2(N), L93) and on the in- 
tensities of mctastable and non-metastable ammonia lines 
(Tkin ~ T mt = 202 ± 15 K in SgrB2(M), Vogel, Genzel & 
Palmer 1987). Note that the introduction of the tiny hot- 
core component does not produce significant variations in 
the predicted continuum emission. 

We also consider the effect of temperature on the 
gas-phase abundance of the modeled species. A vast (~ 
1.7 pc < r < 22.5 pc) portion of the warm envelope shows 
temperatures lower than 100 K, which favour the freeze- 
out of gas- phase water onto dust grains (cf. Williams 1993 
and references therein). The abundance of water (both 
deuterated and non-deuterated) in this region will natu- 
rally be lower than it would be at temperatures higher 
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Table 5. Data set used to model the HDO and H2O distribution in the SgrB2 cloud complex. References are: Jacq 
et aLQ123](J90), Neufeld et al. 120001 (NOP). Gensheimer et al. HOOOI fG96h The set of HDO transitions is displayed in 
Fig. El 



Transition 


Ei/k 


Frequency 


HPBW 


Observed 


Ref. 




(K) 


(GHz) 




towards 




HDO(li,i-U ,o) a 


0.0 


893.6 


11" 


M, N 


this work (CSO) 


HDO(2 1> i-2i, a ) b 


83.6 


241.6 


10" 


N 


J90 (IRAM 30m) 


HDO(3i, 2 -2 2 ,i) b 


156.7 


225.9 


11" 


M, N 


J90 (IRAM 30m) 


HDO(4 2 , 2 -4 2 .3) b 


312.3 


143.7 


17" 


N 


J90 (IRAM 30m) 


o-HfO(l 1 . -lo.i) a 


34.2 


547.7 


3.'9 


M+N 


N00 (SWAS) 


p-H 2 8 0(3 li3 -2 2 ,o) b 


193.9 


203.4 


12" 


N 


G96 (IRAM 30m) 



a absorption 
b emission 



than 100 K, and particularly in the vicinity of the hot 
cores, where the higher T promotes the evaporation of 
water ices from the grains surface. The second major mod- 
ification to the Z95b model involves therefore a differen- 
tiation between a warmer (T > 100 K, hereafter Phase 
I 3 ) and a colder (T < 100 K, hereafter Phase II) region, 
identified by different gas-phase abundances. This is of 
course an oversimplification, however such an approxima- 
tion is sufficient to reproduce the data correctly, within 
the errors, as will be illustrated in § 14.21 As already dis- 
cussed by Neufeld et al. |1997), who introduced a similar 
differentiation in order to reproduce the H 2 0(4 3i2 — 4 2j3 ) 
emission line observed with ISO at 122 /im, this freeze-out 
assumption is required by the observations: a reduction of 
the HDO and H 2 abundance in the outer regions of the 
envelope is necessary to model the HDO and H 2 8 emis- 
sion radiated from the hot-core gas and, at the same time, 
the HDO absorption thought to be produced in the warm 
envelope. In fact, a decrease in the water abundance to- 
wards the outer regions several molecular cloud cores has 
been observed by Snell et al. ( 2000b). The sketch in Fig. 01 
shows the relative sizes of hot core and warm envelope, as 
well as of the Phase I and Phase II regions. 

Finally, as indicated in Tab. all the features in our 
sample but the two o-H 2 8 lines have been observed with 
spatial resolutions ranging between 10" and 17". Since the 
projected distance between the northern and the middle 
core is of ~ 47" (~ 2 pc at a distance of 8.5 kpc), it is more 
appropriate to model the two cores separately. In fact, 
M and N show different chemical compositions, with the 
northern core displaying higher abundances of complex 
molecules (cf. Snyder, Kuan & Miao ll99ll Miao & Snyder 
1997 and references therein; Nummelin et al. 2000). The 
middle core, on the other hand, presents a higher mass 
and luminosity, as inferred from the analysis of the dust 
emission at several wavelengths (cf. Goldsmith et al. [1992 
L93; Gordo n et a l. HOOSI Vogel, Genzel & Palmer 
Dowell et al. 1999). Such dissimilarities seem to reflect the 
youth of N relative to M. However, for our purposes we 

3 Note that Phase I includes, but is not limited to, the hot- 
core component of the cloud (see sketch in Fig. 0J. 



will neglect the chemical differences between the two cores, 
since the absence of molecular species in the middle core 
mostly concerns large complex molecules (see references 
above). Moreover, our simple model shows that the dif- 
ferent continuum temperatures observed towards the two 
cores can be reproduced fairly well by simply assuming dif- 
ferent dust properties for the two cores. L93 suggested, for 
the grain emissivity law in the northern core, a shallower 
slope than that of the middle core, a difference that may 
be explained by different grain shapes in the two cores. 

We proceed by running our radiative transfer model 
separately for the two sources. The results will be appro- 
priately combined, when necessary, to be compared to the 
data (see § I4.2f) . We assume identical physical parame- 
ters (density and temperature profiles, size of the hot core 
and of the warm envelope, H 2 and HDO abundance) 
for both M and N. The only difference between the two 
cores lies in the grain emissivity law, that has a slope of 
1.1 in the northern core, and of 1.4 in the middle one 
(L93). This condition holds only in the inner ~ 1 pc of 
the cloud, which roughly correspond to half the projected 
distance between the two cores. In fact, we must take into 
account the fact that the warm envelopes in which the 
cores are embedded will, at some point, merge to form 
one single "shared" envelope. Hence, for radii larger than 
half the distance between the cores, we assume the dust 
properties to be the same, with a grain emissivity slope of 
1.5 as derived from submm flux densities in a 60" beam 
(Goldsmith et al. nMty . 

4.2. HDO and H2O in the hot cores and in the warm 
envelope 

We aim at reproducing the 143-, 226-, and 241-GHz HDO 
emission features, observed by Jacq et al. |1990| also see 
Tab. |SJ, in order to model the gas-phase abundance of 
this species in the hot-core components of SgrB2(M) and 
N. At the same time, we rely on the 894-GHz ground- 
state absorption feature (Fig. ^) to constrain the HDO 
abundance in the warm envelope. The closest match be- 
tween model and data is obtained when the abundance 
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Fig. 4. Sketch illustrating the relative sizes of the SgrB2 
cloud components. For practical reasons, only one of the 
two cores is represented here. The hot core (labelled HC 
in the figure) has a radius of 0.05 pc, and is tiny compared 
to the size of the warm envelope (WE, r = 22.5 pc). The 
dashed circle identifies the transition radius (r = 1.7 pc) 
between Phase I (T > 100 K) and Phase II (T < 100 K). 
Modeling of the spherically symmetric regions is illus- 
trated in detail in § 14.11 The hot layer (HL) is shown as 
a thin (0.02 pc wide) sheet of gas lying right outside the 
warm envelope (see § I4.3[) . 

of HDO, relative to H2, is set to 1.5 x 10~ 9 in Phase I, 
and 3.5 x 10~ n in Phase II of the molecular cloud, i.e., 
a depletion of a factor of 40 is observed in Phase II. We 
note that this model predicts the 465-GHz line to be in 
weak emission, in agreement with the observations (cf. 
§ |2l ■ The total column density of HDO in a 10" beam, 
iV(HDO), is as high as 2 x 10 16 cm -2 , almost three orders 
of magnitude higher than estimated from the 894-GHz 
absorption line (§ 13.1(1 under the assumption that all the 
absorbing HDO is in the ground state. This inconsistency 
is explained by the radial distribution of the fractional 
populations of the first three levels of HDO (0o,o, lo,i and 
liji in order of increasing energy, Fig. [SJ: the fractional 
population of the ground level exceeds 90% for r <; 3 pc 
only, hence our calculations in S I3.1l severelv underestimate 
the total HDO column density in the innermost portion 
of the cloud, which shows the highest gas density and wa- 
ter abundance and thus contributes the largest percentage 
(> 99%) of the total HDO column density. Incidentally, 
note that our value of TV(HDO) is only about one order 
of magnitude higher than that estimated by Jacq et al. 
(1990) towards SgrB2(N) through the analysis of the hot- 
core transitions, that instead trace the densest gas in the 
complex. 

The abundance of non-deuterated water in the hot 
core can be modeled on the basis of the 203-GHz p- 



HC < > WE Phase I< > Phase II 




10 10 1 10 

Radius (pc) 



Fig. 5. Fractional populations of the three lowest-energy 
levels of the HDO as a function of the radius of the cloud. 



H3; 8 transition observed, in emission, towards SgrB2(N) 
(Gensheimer et al. I1996fl . We are able to reproduce the 
measured intensity of this line if we assume the abundance 
of p-H^ 8 to be 2.2 x 10~ 9 in Phase I. Note that this fea- 
ture shows, towards SgrB2(N), severe blending with a wide 
SO2 line, therefore the intensity indicated by Gensheimer 
et al. should be taken cum grano salis. However, a vari- 
ation of the measured intensity of the 203-GHz line up 
to 50% will not introduce significant changes in our fit 
results. 

Assuming an ortho/para ratio of 3, and a [ 16 0]/[ 18 0] 
ratio of 261±20 (Whiteoak & Gardner 1981), we estimate 
the H2O abundance to be ~ 2.3 x 10~ 6 in this region. 
This result yields a [HDO]/[H 2 0] ratio of 6.5xl0~ 4 , which 
is in order-of-magnitude agreement with the estimate of 
1.8xl0~ 4 by Gensheimer et al. l|1996f) . The discrepancy 
between the [HDO]/[H 2 0] ratio calculated by us, and the 
[D]/[H] ratio measured towards the Galactic Center (cf. 
Lubowich et al. I2000f) . will be discussed in §0 

Knowing the [HDO]/[H 2 0] ratio, the o-H 48 abun- 
dance in Phase II can be determined 4 , [o-H^O] ~ 2.1 x 
10~ 10 . The resulting total peak column density of H 2 
(hot core + warm envelope) is ~ 3.6 x 10 19 cm~ 3 . The 
548-GHz o-H 2 8 absorption expected to be produced 
in the warm envelope is then predicted, for each core, 
with our radiative transfer code. The two spectra are op- 

4 We can reasonably assume that the gas-phase water abun- 
dance at this stage be determined mostly by evaporation from 
dust grains or freeze-out onto dust grains, and that these pro- 
cesses affect both H2O and HDO equally. In other words, we as- 
sume the [HDO]/[H2 0] ratio to be constant throughout Phase 
I and II. However, even considering that fractionation is likely 
to take place in the cold portions of the gas, the [HDO]/[H20] 
ratio in the Galactic Center is expected to assume a maximum 
value of ~ lCT 3 (E. Ber gin, priv. comm.; also cf. Roberts & 
Millar 2000 1 , not very far from our estimate and therefore not 
affecting our conclusions. To make our point even stronger, 
fractionation will be most active in the coldest, outermost re- 
gions of the envelope, which contribute less than 1% to the 
total gas column density. 
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portuncly combined to reproduce the attenuation of the 
SWAS beam, and finally compared to the N00 data. 

A quantitative comparison between the measured in- 
tensities of the whole set of HDO and H3, 8 lines and the 
model results can be found in Tab. Most of the spec- 
tral line data are reproduced within errors of ~ 30%, with 
three exceptions (indicated with a star in Tab.||)|): 

1) The 226-GHz HDO emission line observed towards 
SgrB2(M) is a factor of 2 weaker, and 

2) the 143-GHz HDO line in SgrB2(N) is 4 times stronger 
than predicted by the model. Such discrepancies are 
likely to be due to the intrinsic chemical and physical 
differences between the molecular gas components of 
the two hot cores (cf. Miao & Snvder H997I and refer- 
ences therein) , which have not been taken into account 
in our model. 

3) Most interestingly, the predicted ground-state o-H^O 
absorption at 548 GHz is much shallower than ob- 
served (Fig. [2 N00): the abundance of non-deuterated 
water in the warm envelope is not sufficient to pro- 
duce the absorption feature observed at 548 GHz. 



Table 6. Comparison between the measured intensities 
and opacities of the HDO and H^O lines in our sam- 
ple (Tab. El, and the predictions from the modified Z95b 
model (§ EJ. At this stage, only the hot cores and the 
warm envelope of SgrB2 are included in the model. Note 
that the predicted optical depth of the 548-GHz o-H^O 
line is a factor of 8 smaller than observed (see § 14.21 and 

§123- 



Line 


Source 


Observed T m b 


Model T mb 


(emission) 




(K) 


(K) 


HDO 241 


M 




0.9 




N 


1.2 


0.9 


HDO 226 


M 


0.4 


0.8* 




N 


1.3 


0.8 


HDO 143 


M 




0.1 




N 


0.4 


0.1* 


p-H^ 8 203 


N 


1.5 


1.5 


(absorption) 




Observed r 


Model r 


HDO 894 


M 


0.4 


0.5 




N 


0.8 


0.6 


o-Hf O 548 


M+N 


0.8 


0.1* 



* deviates by more than 50% from observed value 



5 In fact, the column density of o-H3, 8 is of the same order 
of magnitude as that of HDO, but the Einstein B coefficient of 
the o-H2 8 0(li,o — lo,i) transition is about a factor of 2 smaller 
than that of the HDO(li,i — 0o,o) transition, resulting in a 
shallower absorption line. 



4.3. A "hot" issue: H2O in the hot layer 

The result that the observed H^O absorption cannot, ac- 
cording to our model, be produced in the warm envelope 
alone, is a very important one: first of all, our predictions 
confirm, at least from a qualitative point of view, the hy- 
pothesis of C02 that the most important contribution to 
the observed H3, 8 absorption comes from the foreground 
hot gas layer. Secondly, having predicted the contribution 
of the warm envelope to the H^O absorption, it is possible 
to estimate the spatial width and H2O abundance neces- 
sary for the hot layer to produce the observed feature, 
given the observational constraints mentioned in section 
§ 0J In what follows, we will assume the hot layer to be 
located right outside of the warm envelope (see Fig. 0J, 
although its exact location does not matter for our pur- 
poses. 

Because there is no reason to believe it has spherical 
symmetry, the hot layer cannot be modeled directly with 
our radiative transfer code. However, we can use the ra- 
diative transport equation to calculate the total intensity 
emerging from the hot layer, /tot, given a background 
emission Ihc+we as calculated by our model for the com- 
plex made up of hot core and warm envelope (results listed 
in Tab.©: 

/tot = /hc+we • e" THL + /hl • (1 ~ e^), (3) 

where tjjl is the optical depth of the hot layer, and 
/hl • (1 — e _THL ) is the hot layer emission corrected for 
self-absorption. We assume a water abundance, relative 
to H2, of 8 x 10 -4 , based on the abundance of atomic oxy- 
gen not locked in CO (cf. Meyer, Jura & Cardelli 1 1 9981 for 
a determination of the interstellar O abundance; Shaver 
et al. 119831 and Rolleston et al. 120001 for the parametriza- 
tion of the abundance gradient of carbon and oxygen in 
the Galaxy). The spatial width of the hot layer, Aiohl, is 
treated as a free parameter, with 0.01 pc < Awhl < 3 pc, 
and we assume the gas sheet to be extended perpendicu- 
larly to the line of sight (see Fig. HJ|. 

Both dust and gas contribute to thl and /hl - However, 
it is reasonable to assume the contribution of such a thin, 
diffuse layer of dust to be negligible with respect to the 
background envelope, and only the emission and opacity 
of the molecular gas will be taken into account. The frac- 
tional populations of levels l^o and lo,i of o-H3, 8 have 
been calculated, with our radiative transfer code, for a 
thin sheet of hot gas of density tj(H2)hl = 10 3 cm~ 3 
and temperature T gas = 500 K. Both values match the 
lower limits derived by C02 through the modeling of the 
ammonia absorption features observed with ISO at in- 
frared wavelengths. A strict upper limit on the H 2 den- 
sity, n(H 2 )nL Si 5 x 10 3 cm~ 3 , has been determined by 
Hiittemeister et al. |1995J through cm-wavelength obser- 
vations of NH3, SiO and HC3N. These values are further 
supported by measurements of the intensity of the 691- 
GHz CO(6-5) line, performed with the CSO by P. Schilke 
& D.C. Lis (unpublished data). Also, note that the ratio 
of the fractional populations of the o-H^O levels is insen- 
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sitive to temperature changes in the 500-700 K window 
indicated by C02. 

The best fit of the o-H^O column density is deter- 
mined based on the distribution of the \ 2 values given 
by: 

2 (^obs ''"mod) / .\ 

X = f ^ , 4 

where t OD s and r mo d are, respectively, the opacity of the 
observed 548-GHz H5, 8 absorption, and the predicted 
opacity of the feature for each pair of [H^OJhl and Awhl 
values. cr(r obs ) is the error on the observed line opac- 
ity, calculated by adopting a 20% error on the contin- 
uum level measured by SWAS at 548 GHz (N00). The 
minimum \ 2 values yield a o-H3> 8 peak column den- 
sity of ~ 10 14 cm~ 2 , hence a H2O column density of 
~ 3.5 x 10 16 cm~ 2 , corresponding to Awhl ~ 0.02 pc. 
We consider the value of Awhl to be a lower limit, but 
expect it to be correct within a factor of a few (~ 2-3). 
If we assume a [HDO]/[H 2 0] ratio equal to the [D]/[H] 
ratio measured towards the Galactic Center (1.7 X 10~ 6 , 
Lubowich et al. 2000) as expected for gas at this tem- 
perature, then the HDO peak column density in the hot 
layer, A(HDO) H l ~ 5.2 x 10 10 cm~ 2 , will contribute with 
0.3% only to the absorption produced by the HDO com- 
ponent located in the warm envelope, in agreement with 
our hypothesis (§EJ. 

An overview on the radial density and temperature 
profiles adopted to model the three regions of the cloud 
complex is given in Tab.0 Each profile, P(r), is described 
by the equation 

P(r) = a + a 1 ■ . (5) 

Also, Tab. [S] summarizes the peak column densities of H 2 , 
HDO and H 2 resulting from our model. The H 2 abun- 
dance and the [HDO]/[H 2 0] ratio are also stated. 

5. Discussion 

The unambiguous identification of the hot layer as the gas 
component responsible for the 65-kms -1 water absorp- 
tion towards SgrB2 is definitely the main result of this 
work. This confirms the general picture proposed by C02. 
However, our model provides, for the hot layer, a value of 
the water column density, iV(H 2 0)HL, which is a factor 
of ~ 10 lower than that inferred by C02, due to the fact 
that we assumed a lower H 2 density. In fact, previous es- 
timates of the column density of the absorbing water (cf. 
Z95a, N00, and ? !3.2l of this paper) are in very good agree- 
ment with our result of iV(H 2 0) H L ~ 3.5 x 10 16 cm" 2 : 
the assumption employed in these estimates that only the 
ground state of H 2 be populated is correct for the hot 
layer, although it does not apply to the warm envelope. 
However, the column density of water actually responsible 
for the observed absorption is negligible, with respect to 
the total column density estimated over the whole enve- 
lope (see Tab. IBJ. It is important to keep in mind that 



column densities calculations from ground-state absorp- 
tion features, if based on the "classical" assumptions men- 
tioned in § 13.11 can be highly misleading in regions as 
complex as SgrB2, and are likely to yield underestimates 
of several orders of magnitude (a factor of ~ 1000 in our 
case) . 

As previously mentioned, several heating mechanism 
have been proposed to explain the existence of such a hot, 
relatively thin sheet of molecular gas. The origin of the 
whole cloud complex has been proposed to be linked to 
large-scale cloud-cloud collisions (Hasegawa et al. D-994), 
so shock-induced heating is definitely possible, although it 
remains to be explained why the velocity of the hot layer 
is identical to that of the warm envelope. Martin-Pintado 
et al. il2lT(l()i argue rather in favour of an X-ray-driven 
chemistry, based on the observed spatial correlation of 
the Fe 6.4-keV emission line with the SiO(l-O) emission at 
43.4 GHz. In a recent paper by Goicoechea & Cernicharo 
(2002), a [OH]/[H 2 0] abundance ratio of 0.1-1 is estimated 
for the hot layer, which, according to the authors, points 
towards the presence of a strong UV field illuminating the 
outer shells of the cloud. The high abundance of H 2 es- 
timated by us does not allow to discriminate among the 
proposed mechanisms, and it is anyway likely that all of 
them contribute, to some extent, to the anomalous heating 
of this region. However, our radiative transfer calculations 
allow us to separate the chemistry driving the water abun- 
dance in the hot layer from that taking place in the warm 
envelope, thus helping to set more solid constraints on the 
physical characteristics of the hot diffuse gas, such as, for 
example, its spatial extent. 

We would also like to stress the high deuterium frac- 
tionation ([HDO]/[H 2 0] ~ 6.4 x 10~ 4 ) inferred, for the 
SgrB2 hot cores, by our radiative transfer model. Our 
value is in agreement (about a factor of 3 higher) with that 
estimated by Gensheimer et al. (1996) for SgrB2(N), and 
it is almost 400 times larger than the elemental [D]/[H] ra- 
tio measured towards the Galactic Center 6 (Lubowich et 
al. [?000 ) . An enhancement of the deuterium fractionation 
in hot-core-type sources has been observed for a variety 
of chemical species ([CHaODj/pHjOH], Mauersberger et 
aLUHHEl [NH 2 D]/[N H 3 ], Jacq et al. lTMH [DCN]/[ HCN], 
Schilke et al. HM?1 Hatchell, Millar & Rodgers UW® 
[DCO+]/[HCO+], Jacq e t al. HTO [HDS]/[H 2 S], Hatchell, 
Ro berts fe Millar ITWfll [CH 2 DCN]/[CH 3 C N], G e rin et 
al. and [ HDO]/ [H 2 Q], Jacq et a l. ITMH 17333 

Gensheimer et al. HHHU Pardo et al. EDDU, and is gener- 
ally attributed to evaporation of deuterated species from 
grain mantles due to the formation of an embedded heat- 
ing source (see Walmsley et al. I1987|) . It must be kept 
in mind that the quoted estimates of the [D]/[H] ratio 
in the Galactic Center are affected by large uncertainties 
(about 1 order of magnitude). However, we can assume 



The ratio of the abundance of a deuterated species with 
respect to its non-deuterated counterpart, [XD]/[XH], is ex- 
pected to scale with the [D]/[H] ratio in the gas (Gerin & Roueff 
119991 E. Bergin, priv. comm.). 
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Table 7. Description of the radial profiles of density and temperature used to model the hot cores and warm envelope 
(HC and WE respectively, § 14.2(1 . and the hot layer (HL, see § 14.3(1 of the SgrB2 cloud complex. The profiles are 
described by the coefficients a , a\, r and a according to the law: P(r) = a + a\- (r/r ) a . The radii are (cf. Fig.QJ: 
r < 0.05 pc for the hot core and 0.05 pc < r < 22.5 pc for the warm envelope. The hot layer is, in our model, 0.02 pc 
thick. 



n(&2) Tgas Tdust 

ao ai ro ct ao a\ ro ct ao ai ro a 

Region (cm" 3 ) (cm" 3 ) (pc) (K) (K) (pc) (K) (K) (pc) 

HC 3.4 x 10' 1 ~ 200 - - 200 

WE 2.2 x 10 3 5.5 x 10 4 1.25 - 2 40 1 -0.5 40 1 -0.5 

HL 10 3 500 



Table 8. Summary of the peak column densities of H2, HDO and H2O derived in § 14.21 and 14.31 The values relative 
to hot cores and warm envelope (HC+WE) are based on the radiative transfer code illustrated in Z95b and modified 
as in § 14.11 The determination of the hot layer values is described in S 14.31 The [HDO]/[H2 0] ratio for each region is 
also indicated. The value of [HDOJ/prbO] in the hot layer is assumed equal to the [D]/[H] ratio in the Galactic Center 
(cf. Lubowich et al. lMirjji . 



Region 


iV(H 2 ) 
(cm" 2 ) 


iV(HDO) 

(cm" 2 ) 


7V(H 2 0) 
(cm" 2 ) 


[HDO]/[H 2 0] 


HC+WE 
HL 


1.9 x 10 2b 
6 x 10 19 


2.3 x 10 10 
5.2 x 10 10 


3.6 x 10 iH 
3.5 x 10 16 


6.4 x 10-* 
1.7 x 1Q- 6W 



W assumed (cf. Lubowich et al. 2000 ) 



our result to be in agreement with the general finding 
that the abundance ratio of deuterated species to their 
non-deuterated counterparts is, in hot-core-type sources 
(cf. references above), enhanced by a factor of a few 10 2 
with respect to the elemental [D]/[H] ratio. Thus, the value 
of [HDO]/[H 2 0] derived in §Q yields a [D]/[H] ratio of 
a few 10" 6 , supporting the hypothesis of a lower [D]/[H] 
ratio in The Galactic Center region. 
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